Production cross sections from 82 Se fragmentation as indications 
of shell effects in neutron-rich isotopes close to the drip-line 
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Production cross sections for neutron- rich nuclei from the fragmentation of a 82 Se beam at 
139 MeV/u were measured. The longitudinal momentum distributions of 126 neutron-rich iso- 
topes of elements 11 < Z < 32 were scanned using an experimental approach of varying the target 
thickness. Production cross sections with beryllium and tungsten targets were determined for a 
large number of nuclei including several isotopes first observed in this work. These are the most 
neutron-rich nuclides of the elements 22 < Z < 25 ( 64 Ti, 67 V, 69 Cr, 72 Mn). One event was regis- 
tered consistent with 70 Cr, and another one with 75 Fe. The production cross sections are correlated 
with Q g systematics to reveal trends in the data. The results presented here confirm our previous 
result from a similar measurement using a 76 Ge beam, and can be explained with a shell model that 
predicts a subshell closure at N — 34 around Z — 20. This is demonstrated by systematic trends 
and calculations with the Abrasion- Ablation model that are sensitive to separation energies. 

PACS numbers: 25.70.Mn, 27.40.+Z, 27.50.+e, 21.60.Cs 



I. INTRODUCTION 
A. Discovery of new nuclei 

The discovery of new nuclei in the proximity of the neu- 
tron dripline provides a stringent test for nuclear mass 
models, and hence for the understanding of both the nu- 
clear force and the creation of elements. Another impor- 
tant aspect of such measurements is that once neutron- 
rich nuclei are observed and their cross sections for for- 
mation are understood, investigations to study the nu- 
clei themselves, such as with decay spectroscopy, can be 
planned. Therefore, obtaining production rates for the 
most exotic nuclei continues to be an important part of 
the experimental program at existing and future rare- 
isotope facilities. 

A number of production mechanisms have been used 
to produce neutron-rich isotopes for 20 < Z < 28 [l| 
but, in the last few years, two reaction mechanisms were 
the most effective at producing nuclei in this region: 

• projectile fragmentation - an experiment with a 
76 Ge (132 MeV/u) beam produced 15 new isotopes 
of 17 < Z < 25 0, 

• in-flight fission with light targets (Abrasion- 
Fission) - an experiment with a 238 U beam Q pro- 
duced a large number of isotopes of 25 < Z < 48 
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using a Be-target, and several new isotopes with 
46 < Z < 56 by Coulomb fission on a heavy target. 

Progress in the production of neutron- rich isotopes was 
made possible by the increase of primary beam intensi- 
ties, new beam development at the National Supercon- 
ducting Cyclotron Laboratory (NSCL) at Michigan State 
University and advances in experimental techniques Q. 
Indeed, recent measurements at the NSCL P, U@ have 
demonstrated that the fragmentation of 48 Ca and 76 Ge 
beams can be used to produce new isotopes in the prox- 
imity of the neutron dripline. Continuing this work, we 
report here the next step with a newly developed 82 Se 
beam towards the fundamental goal of defining the ab- 
solute mass limit for chemical elements in the region of 
calcium. In the present measurement, four neutron-rich 
isotopes with 42 < TV < 47 were identified for the first 
time (see FigJT]), one event was registered consistent with 
70 Cr46, and another one with 75 Fe4g. 



B. Evidence for global structure changes 

One of the first indications of significant changes in 
the structure of neutron rich nuclei was the discovery of 
enhanced nuclear binding of heavy sodium isotopes [|| . 
This is now understood to result from significant contri- 
butions of fp shell intruder orbitals to the ground-state 
configurations of these isotopes 0, [Io|. Low- lying 2 + 
states and quadrupole collectivity have been reported 
in neutron-rich even-even Ne and Mg isotopes around 
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FIG. 1. (Color online) 
The region of the nu- 
clear chart investigated 
in the present work. 
The solid line shows the 
limit of bound nuclei 
from the KTUY mass 
model Q. The new iso- 
topes observed for the 
first time in the present 
work are marked by red 
squares. 



N = 20, see for example Refs. [TT| - [l5| . This region 
around 31 Na, where the neutron fp shell contributes sig- 
nificantly to the ground-state structure, is now known as 
the "Island of Inversion" . Similarly, there is mounting ev- 
idence for an onset of deformation around neutron num- 
ber N ~ 40 in Fe and Cr nuclei. In even-even Fe and Cr 
nuclei, for example, this evidence is based on the energies 
of low- lying state s fl6l - [2rjj , transition strengths [2l| , de- 
formation length |l9|. and higher-spin level schemes f22j . 
Neutron <7g/ 2 and configurations from above the 

N = 40 shell gap arc proposed to descend and domi- 
nate the low-lying configurations similar to those in the 
N = 20 Island of Inversion [H, [H[. 

Recently, it was shown that the beta-decay half-lives of 
the neutron-rich Ca isotopes [25j compare favorably with 
the results of shell-model calculations performed in the 
full pf model space using the GXPF1 effective interac- 
tion [26[ . The systematic trend of these half- lives is con- 
sistent with the presence of a subshell gap at N = 32 as 
predicted by this interaction and confirmed by a variety 
of experiments. This interaction also predicts an increase 
of the excitation energy of the first excited state E x (2f) 
at 54 Ca relative to that obtained with the KB3G [27fl 
interaction, suggesting the appearance of the N = 34 
shell gap in Ca isotopes. Both interactions predict similar 
structure for light stable nuclei, but give rather different 
predictions for several cases of neutron-rich nuclei. 

Recent measurements at RIKEN of the E X {2 J [) in 
54 Ca 0| found that the experimental value is 0.5 
MeV smaller than the GXPF1B prediction, where 
the GXPF1B [29| Hamiltonian was created from the 
GXPF1A Hamiltonian by changing five T = 1 matrix 
elements and the single-particle energies which involved 



lpi/2- A similar trend had already been pointed out by 
Mantica et al. [25| where they deduce that the effective 
energy gap between the adjacent neutron single-particle 
orbitals / 5 / 2 and pi/ 2 is overestimated by the GXPFI 
and GXPFIA effective interactions. Based on this, the 
GXPFIB interaction has been modified to correct this 
0.5 MeV shift and is referred to as GXPFIB5 here. 



Other evidence supports the modified form of the 
GXPFIB interaction. The original GXPFIB interaction 
predicts a ln-unbound 56 K (S ln = -0.03 MeV); how- 
ever, this isotope was shown to be bound by observation 
in our previous experiment with the 76 Ge beam. The 
shift in the interaction makes the isotopes with valence 
neutrons in the / 5 / 2 orbital around Z = 20 more bound, 
such that the modified interaction GXP1FB5 predicts a 
bound 56 K with S ln = 0.41 MeV. 



In our previous cross section measurements in the re- 
gion around 62 Ti ( 76 Ge primary beam) [2j we observed 
a systematic variation of the production cross sections 
that might point to nuclear structure effects, such as an 
onset of collectivity, that are not included in global mass 
models that were used to construct the basis of the sys- 
tematics. The present work, since it is based on isotope 
production from a different primary beam, covering the 
same region of the nuclear chart, provides an independent 
check of this interpretation. 
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TABLE I. Experimental settings 
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II. EXPERIMENT 
A. Setup 

A newly developed 139 MeV/u 82 Se beam with an in- 
tensity of 35 pnA, accelerated by the coupled cyclotrons 
at the NSCL, was fragmented in a series of beryllium 
targets and a tungsten target, each placed at the object 
position of the A1900 fragment separator [3(j. In this 
work we used an identical configuration to our previous 
experiment with a 76 Ge beam [1[, where the combination 
of the A1900 fragment separator with the S800 analysis 
beam line [3l| formed a two-stage separator system, that 
allowed a high degree of rejection of unwanted reaction 
products. At the end of the S800 analysis beam line, 
the particles of interest were stopped in a telescope of 
eight silicon PIN diodes (50x50 mm 2 each) with a total 
thickness of 8.0 mm. A 50 mm thick plastic scintillator 
positioned behind the Si-telescope served as a veto de- 
tector against reactions in the Si-telescope and provided 
a measurement of the residual energy of lighter ions that 
were not stopped in the Si-telescope. A position sensitive 
parallel plate avalanche counter (PPAC) was located in 
front of the Si-telescope. All experimental details and a 
sketch of the experimental setup can be found in Ref. [l[ . 
In this paper, we describe the details of our experimental 
approach and discuss the results. 



B. Experimental runs 

The present experiment consisted of four segments that 
are summarized in Table HI Except for the last segment, 
the present experimental program is similar to the previ- 



ous Ge experiment During all runs, the magnetic 
rigidity of the last two dipoles of the analysis line was 
kept constant at a value of 4.2867 Tm while the produc- 
tion target thickness was varied to map the fragment mo- 
mentum distributions. This approach greatly simplifies 
the particle identification during the scans of the parallel 
momentum distributions. 

The momentum acceptance of the A1900 fragment sep- 
arator was restricted to Ap/p = 0.1% (first four runs with 
thin targets), and to Ap/p = 0.2% (other targets) for 
the measurement of differential momentum distributions 
in the first part of the experiment. The use of different 
beryllium target thicknesses (9.7, 68, 138, 230, 314, 413, 
513 mg/cm 2 ) allowed coverage of the fragment momen- 
tum distributions necessary to extract production cross 
sections and also resulted in more isotopes in the particle 
identification spectrum. 

For the second part of the experiment, a Kapton wedge 
with a thickness of 20.0 mg/cm 2 was used at the disper- 
sive image of the A1900 with a 10 mm aperture in the fo- 
cal plane to reject less exotic fragments while the separa- 
tor was set for 67 Fe and 78 Zn ions. The goal of this setting 
was to confirm the particle identification by isomer tag- 
ging as described in Ref. H3 with 67m Fe (E y = 367 keV, 
T 1/2 = 43 /its) and 78m Zn [E 1 = 730,890,908 keV, 
Ti /2 = 0.32 /is). 

In the third part of the experiment, dedicated to the 
search for new isotopes, five settings were used to cover 
the most neutron-rich isotopes with 20 < Z < 27, as it 
was impossible to find a single target thickness and mag- 
netic rigidity to observe all of the fragments of interest. 
Each setting was characterized by a fragment for which 
the separator was optimized. A search for the most ex- 
otic nuclei in each setting was carried out with Be and 
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A-3Z 

FIG. 2. (Color online) Particle identification plot showing the 
measured atomic number, Z, versus the calculated function 
A — 3Z for the nuclei observed in production runs of this work. 
See text for details. The limit of previously observed nuclei is 
shown by the solid red line, and the location of 64 Ti is marked. 



W targets. The settings were centered on 60 Ca, 68 V and 
74,75p e respectively, based on LISE ++ [33| calculations 
using the paramctcrizations of the momentum distribu- 
tions obtained in the first part of the experiment (see Sec- 
tion |IV3}. The momentum acceptance of the A1900 was 
set to the maximum of Ap/p = 5.0% for these production 
runs. It should be noted that the momentum acceptance 
of the S800 beamline is about 4% according to LISE ++ 
Monte Carlo simulations using a new extended configu- 
ration with 5 order optics. This calculated acceptance 
has been used for the cross section analysis using the 
method described below. 

The fourth part of the experiment was devoted to two 
short runs measuring the yields of more stable isotopes 
by centering on 45,48 Ca. 



III. ANALYSIS OF EXPERIMENTAL DATA 

The result of our approach of keeping the last dipoles 
magnetic rigidities constant while varying the target 
thickness — as was done in the previous experiment — 
can be seen in Fig. [21 which shows the total distribu- 
tion of fully-stripped reaction products observed in the 
production runs of this work. The range of fragments 
is shown as a function of the measured atomic number, 
Z, versus the quantity A — 3Z deduced from measured 
values, where A is the mass number. The identification 
of the individual isotopes in Fig. [5] was confirmed via 
isomer tagging using the known isomeric decays in 67 Fe 
and Zn. The standard deviations of ionic charge (<?) 
and elemental (Z) spectra were found to be similar to 
those in the previous experiment, therefore the proba- 
bilities of one event being misidentified as a neighboring 
charge state or element as before [![. The details of the 
calculation of the particle identification are given in the 
appendix to the previous work [l|. 

The mass spectra for the isotopic chains from scandium 
to iron measured during the production runs are shown 
in Fig. [3] Only nuclei that stopped in the Si telescope 
are included in this analysis. The observed fragments 
include several new isotopes that are the most neutron- 
rich nuclides yet observed of elements 22 < Z < 25 ( 64 Ti, 
67 V, 69 Cr, 72 Mn). One event was found to be consistent 
with 70 Cr, and another one with 75 Fe. The new neutron- 
rich nuclei observed in this work lie to the right of the 
solid line in Fig. [5] and to the right of the vertical dashed 
lines in Fig. [3] 

IV. RESULTS AND DISCUSSION 
A. Parallel momentum distributions 

The prediction of the momentum distributions of 
residues is important when searching for new isotopes 
in order to set the fragment separator at the maximum 
production rate. Also, the accurate prediction of the 
momentum distributions allows for a precise estimate 
of the transmission and efficient rejection of strong con- 
taminants. In this experiment the "target scanning" ap- 
proach (3~H ]. developed in the previous experiment, was 
used to obtain parameters for the neutron-rich isotope 
momentum distribution models such as [35l [36j . This 
method is particularly well suited to survey neutron-rich 
nuclei since the less exotic nuclei are produced with the 
highest yields and their momentum distributions can be 
readily measured with thin targets. 

The data analysis of this approach has been improved, 
and a detailed description is in preparation [33]. Im- 
portant improvements include: first, the most probable 
velocity for a fragment is not that at the center of the 
target when the yield is sharply rising or falling with mo- 
mentum, and second, asymmetric Gaussian distributions 
have been used with asymmetry coefficients taken from 
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FIG. 3. (Color on- 
line) Mass spectra 
of the elements 
21 < Z < 26 that 
were stopped in the 
Si-telescope during 
the production runs. 
The limits of previ- 
ously observed nuclei 
are shown by the 
vertical dashed lines. 
Standard deviations 
produced with the 
Gaussian function 
at constant width 
(dashed curves) are 
given in the panel 
figures for each ele- 
ment. 



Mass number 



the convolution model |38( implemented in the LISE ++ 
code [33[. Note that, at the bombarding energy used in 
these experiments, the shape of the fragment momentum 
distribution is asymmetric with a low- energy exponential 
tail stemming from dissipative processes [381 ] . 

Seven targets were used to measure the momentum dis- 
tributions (see Table fl} . The yield of one example frag- 
ment, 66 Fe, is shown in Fig.0]as a function of the ratio of 
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FIG. 4. (Color online) Differential cross sections of 66 Fe frag- 
ments obtained by varying the target thickness with the anal- 
ysis line set at one fixed magnetic rigidity. The dashed line 
represents the fitted asymmetric Gaussian function. Blue hor- 
izontal errors with diamonds in center correspond to the ve- 
locity difference caused by production at the beginning or the 
end of target, whereas the red circles show the position of the 
most probable velocity based on the momentum distribution 
parametrization. 



fragment and beam velocities. This figure illustrates the 
impact of the new data analysis where the most probable 
values taken are shown by the circles and the average val- 
ues by the diamonds. Momentum distributions for 126 
isotopes were derived (indicated by the colored boxes in 
Fig. [T|) and integrated to deduce the production cross 
sections. 

A survey of all of the fitted results showed that frag- 
ments in the heavy mass region were produced simi- 
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FIG. 5. (Color online) Widths of parallel momentum com- 
ponent as a function of mass number of fragments produced 
in the reaction 82 Se beams with beryllium targets. Small 
diamonds denote calculations by the convolution model [38|] 
with default settings for separation energy (E s ) option #1 in 
LISE ++ . Solid green and dot-dashed black lines represent the 



best fit to the data for the Goldhaber 
models, respectively. 



[35] and Morrissey [36 
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lar to our previous measurements [341 ] with significantly 
higher velocities than the momentum distribution mod- 
els [3(| HsJ predict. The difference is most likely due to 
the fact that the models were developed for fragments 
close to stability, where the energy required to remove 
each nucleon was set to 8 MeV, while the actual nucleon 
binding energy for the neutron rich isotopes under in- 
vestigation is lower. An analysis with asymmetric dis- 
tributions to reproduce the mean velocity of fragments 
has shown that the neutron-rich separation energy pa- 
rameter in the model [36| for the nuclei observed in the 
present work in the region Ap/2 < Ap < Ap can be rep- 
resented by a linear decrease with the number of removed 
nucleons: 

E s = 8 - 9.2AA/A P (1) 

where A A = Ap — Ap, Ap is the projectile mass number, 
and Ap is the fragment mass number. In the fourth part 
of the experiment, where stable isotopes were measured 
(see Table H|, no deviations from the default parameters 
of the model for the velocities were observed. 

The width obtained for each fragment's parallel mo- 
mentum distribution is presented in Fig. [3] for fragments 
produced from the interaction of 82 Se (139 MeV/u) with 
9 Be targets. The predictions with best fits of reduced 
widths cr (G) = 87 MeV/c for the Goldhaber [H and 
<7o(M) = 75 MeV/c for the Morrissey [36| models are 
presented in this figure. 

B. Production cross sections 

The inclusive production cross sections for the ob- 
served fragments were calculated by correcting the mea- 
sured yields for the finite momentum and angular ac- 
ceptances of the separator system. A total of 126 cross 
sections with beryllium targets were obtained from the 
Gaussian fits to the longitudinal momentum distribu- 
tions; these nuclei are indicated by stars in Fig. [5] The 
cross sections for all of the remaining fragments with in- 
completely measured longitudinal momentum distribu- 
tions were obtained with estimated transmission correc- 
tions as in our previous work [l[ . The parameters for the 
transmission corrections were assumed to be smoothly 
varying with A and Z . 

The cross sections obtained for all the fragments ob- 
served in this experiment are shown in Fig. [6J along 
with the predictions of the recent EPAX 3 parameteriza- 
tion |4(|. For those isotopes that relied on transmission 
calculations, the weighted mean of all measured yields 
was used to obtain the cross section (shown by solid dia- 
monds in Fig. The uncertainties in these cases include 
the statistical, the systematic and the transmission cor- 
rection uncertainties. For more details see ref. HH . As 
can be seen in Fig. [6j the cross sections are in good agree- 
ment with those produced by integrating the measured 
longitudinal momentum distributions in the cases where 
there is an overlap. 



It is important to note that the predictions of the re- 
cent EPAX 3 parameterization for reactions with beryl- 
lium, shown by the solid lines in Fig. [51 reproduces the 
measured cross sections for isotopes much better than 
the previous EPAX 2.15 predictions j4lj . 

C. Q 9 systematics 

The production cross sections for the most neutron- 
rich projectile fragments have been previously shown to 
have an exponential dependence on Q g , where Q g is de- 
fined as the difference in mass-excess between the beam 
particle and the observed fragment @, @- To test this 
behavior, the cross sections for each isotopic chain were 
fitted with the simple expression: 

a(Z, A) = k(Z) exp (Q g (Z, A))/T{Z)), (2) 

where T represents an inverse slope parameter, and k is 
a normalization. 

Most of the data from the reactions of 82 Se on Be tar- 
gets in this experiment could be fitted by two inverse 
slopes with a floating connection point. The trends of 
the general increase in T for all of the heavy isotopes of 
elements Z = 19, 20, and 21 observed with a 76 Ge beam 
in our previous experiment is reproduced here with the 
82 Se beam. The cases of Z = 16, 19, 20, and 22 are il- 
lustrated in Fig. [7] showing the measured cross sections 
versus Q g calculated using the masses deduced from the 
shell model with the GXPF1B5 interaction [2<|. As in the 
previous experiment, the heaviest isotopes of elements in 
the vicinity of Z = 20 appear to deviate from an expo- 
nential dependence. The change of slope is most dra- 
matic at Z — 20. At Z — 22, the trends show little 
change in slope for heaviest isotopes. As was pointed 
out in our previous study 0], a possible explanation for 
the exponential slope reduction at larger masses for these 
elements is a change in binding energy relative to predic- 
tions. It should be noted that the systematic variation 
of the production cross sections of the calcium isotopes 
as a function of Q g was checked in our previous work 
with several other well-known mass models and essen- 
tially the same behavior was observed (for details, see 
Fig. 10 in 

A hint at the origin of the cross section systematics 
may be seen in the binding energy trends demonstrated 
in Fig. [8j where the dependence of Szn is shown as a func- 
tion neutron number in the neutron-rich region for sulfur, 
potassium, calcium and titanium based on masses from 
models and experimentally measured values. There is no 
abrupt behavior for potassium and calcium in the KTYU 
model or the AME2012 mass table [Hj]. The slopes of 
the trends for both elements in the region 34 < N < 40 
do not change significantly. On the other hand for two- 
neutron separation energy lines calculated with the full 
p/-shell phenomenological GXPF1B5 (2(| interactions, 
changes in the slope are observed between N = 35 and 
N = 36 for potassium and calcium isotopes near the 
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FIG. 6. (Color online) Inclusive production cross sections for fragments from the reaction of Se with beryllium and tungsten 
targets shown as a function of mass number. The cross sections with the beryllium targets derived by momentum distribution 
integration are shown by stars, those normalized with LISE ++ transmission calculations are indicated by solid diamonds The 
cross sections obtained with the tungsten target were normalized with LISE ++ transmission calculations. The red solid lines 
show the predictions of the EPAX 3 systematics [iTj] for beryllium (see text). The two magenta dashed lines separate nuclei that 
require neutron pickup in the production mechanism. 



same neutron number where the cross section systemat- 
ics change. Because 55 Ca (iV=35) is predicted by the 
shell-model with the GXPF1B5 effective interaction [29| 
to have a low one- neutron separation energy of 0.75 McV, 
we might expect that the change of slope in Fig. [7] would 
correspond to N = 35 instead of N = 36. Thus, it is pos- 
sible that the trends seen in our cross section data reflect 
the N = 34 subshell closure predicted by the GXPF1B5 
shell model [2{|. Similar breaks in the slopes are seen in 
the data at N = 28, but they are not as dramatic as at 
N = 34. 



D. Global trends of cross sections 

A way to visualize the possible effect of the N = 34 
subshell closure is to plot the production cross sections 
versus atomic number. Fig. [5] shows the data connected 
by lines of constant N — 2Z, which represent changes be- 
tween nuclei different by a triton, and the label at each 
point is the neutron number (N). This ensemble of lines 
exhibits a large dip at the shell closure at Z = 20 in the 
region of 54-56 Ca highlighted by the red-dashed quad- 
rangle. The same dip can be observed in a similar figure 
generated from the cross sections measured in the 76 Ge 
measurements 

A reason that these trends may be visible in lines of 
TV — 2Z in Fig. [HI is that such nuclei have approximately 
the same neutron separation energy and the drip-line lies 
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FIG. 7. (Color online) Measured cross sections versus -Q g for 
the production of neutron-rich calcium (top plot) and potas- 
sium (middle plot) isotopes from reactions of 82 Se and 76 Ge 
with a beryllium target. Production cross section of sulfur and 
titanium isotopes from the reaction of 82 Se with a beryllium 
target are shown in bottom panel. See text for explanation 
of Q g and the lines. Neutron numbers of isotopes are shown 
in figure. A horizontal dash green line in top plot shows the 
cross section limit reached in the experiment with the 82 Se 
beam. 
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FIG. 8. (Color online) The two-neutron separation energy 
S2n deduced from mass values as function of neutron number 
for calcium (top), potassium (middle), sulfur and titanium 
isotopes (bottom). Those from experimental mass values [i^] 
are shown by diamonds and from AME2012 [i^] by crossed 
circles. Results based on the full p/-shell phenomenological 
GXPF1B5 [H] interaction and the KTYU mass model @] 
are shown by solid squares and hollow triangles, respectively. 
Arrows show regions of isotopes whose measured cross section 
are shown in Fig. [7] 



close to an N — 2Z line for 16 < Z < 23. Further, each 
line connects nuclei with either an odd or even number of 
neutrons; hence, the large odd-even difference in nuclear 
binding due to pairing is not present along each line. For 
a constant separation energy we expect the cross section 
to fall smoothly with Z in Fig. HI At Z = 20, N = 34 the 



(sub)shell closures correspond to a lower S\ n and 5^ and 
might be responsible for the dip in the trends, highlighted 
by the red box. 

To illustrate this, two-neutron separation energy S2 n 
versus neutron number (N) of elements 12 < Z < 22 cal- 
culated with the full p/-shell GXPF1B5 interactions ^ 
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FIG. 9. (Color online) Production cross section versus atomic 
number (Z) for fragments from reaction of 82 Se with beryl- 
lium targets. Lines are connected according to constant 
N — 2Z, while labels represent the neutron number. Re- 
actions resulting in neutron pick-up are omitted. The red 
dashed quadrangle is explained in the text. 



FIG. 10. (Color online) Two-neutron separation energy S2n 
versus neutron number (N) for elements 12 < Z < 22. Values 
are calculated using results from the GXPF1B5 [2!| model. 
Labels in the lines show atomic numbers of nuclei. The red 
dashed quadrangle is explained in the text. 



E. Abrasion-Ablation model 



is shown in Fig. [TO] The label at each point is the atomic 
number (Z) and the red-dashed quadrangle encompasses 
the same region as the one in Fig. [9] for isotopes with 
34 < N < 37 and 18 < Z < 21. The N = 34 sub- 
shell closure results in a lowering of the S^n values for 
Z = 20 as seen in the figure, corresponding to the same 
nuclei that have relatively low production cross sections 
compared to the N — 2Z trends in Fig. [9l 

Based on the observations discussed above, it seems 
possible that the production cross section systcmatics 
provide a hint of a change of the nuclear mass surface. 
The effect is most pronounced close to the drip- line. Plot- 
ting the cross section of elements against Q g can exhibit 
sudden changes in slope that are correlated with regions 
of changes in the nuclear structure, such as (sub)shcll 
closures. 



In order to test the relationship between cross sections 
and separation energies, the production cross sections ob- 
tained in this experiment were compared with those from 
calculations with the Abrasion- Ablation model [44| im- 
plemented in the LISE ++ code [33j]. Results for neutron- 
rich isotopes of potassium and calcium using different 
mass models are shown in Fig. [TT] Predictions of the 
EPAX 3 systematics are shown on the plots for compar- 
ison. An excitation energy of 15 MeV per abraded 
nucleon was deduced from an experimental data fit 
of the Abrasion-Ablation (AA) model with AME2012 
masses [43[ and masses deduced from the shell model 
using the GXPF1B interactions [26j], whereas a value 
18 MeV has been obtained with KTYU masses Q- 
The LISE ++ LDM1 paramaterization was used to ex- 
trapolate masses of very neutron-rich nuclei absent in 
AME2012 and GXPF1B mass predictions. The AA 
model with GXPF1B masses significantly underestimates 
cross sections for isotopes with iV > 35 such as 54 > 55 K 



10 




25 30 35 25 30 35 40 

Neutron Number (N) 



FIG. 11. (Color online) Top: The cross section versus neutron number for the production of neutron-rich potassium (left) 
and calcium (right) isotopes from the reaction of 82 Se on beryllium targets. The measured values are shown by semi-solid 
rectangles. The results from the Abrasion- Ablation calculations using the masses obtained via both the GXPF1B and GXPF1B5 
interactions [2tJ are shown by open diamonds with dashed lines and solid lines correspondingly. For comparison, the values 
obtained from the EPAX 3 systematics [13] (small solid rectangles with dashed-dot-dot lines) are also included. 
Bottom: Logarithm ratios of calculated and experimental cross sections shown for clearer presentation differences between 
experimental and calculated values. This presentation includes results from the KTYU model [Tfl (solid triangles with dot lines) 
and the AME2012 table [43j] (semi-solid circles with dash-dot lines). 



and 55-57 Ca, whereas fair agreement is observed us- 
ing GXPF1B5 masses. Decreasing the effective energy 
gap between adjacent neutron single-particle orbitals f§/2 
and P1/2 in GXPF1B makes isotopes with 36 < N < 40 
around calcium more particle bound, therefore the result 
of the 54 Ca E x (2f) measurement [28( is in good agree- 
ment with the measured cross sections and whose calcu- 
lated using the AA model with GXPF1B5 masses. 



On the other hand the mass models that don't pre- 
dict the slope changes at N = 36 in the S^n figure 
(see Fig. [5]) overestimate the cross sections for neutron- 
rich potassium and calcium isotopes (see AA calculations 
with masses from the KTYU model |7fl in Fig. [TT]) . 



F. Estimation of Ca production cross section 

Based on experimental Q g -systematics (FigJT]) and 
Abrasion- Ablation model calculations (Fig. ITTj) using 
GXPF1B5 masses, it is possible to estimate the produc- 
tion cross sections for the next unobserved calcium iso- 
topes using a 82 Se beam (See Tabic HI]) . According to this 
extrapolation, to observe one event of 59 Ca, the beam in- 
tensity has to be increased by a factor of 2 compared to 
this experiment. For the case of 60 Ca, at least an or- 
der magnitude higher beam intensity is needed. It is 
important to note, that the EPAX 3 systematics pre- 
dict a factor of 20 less production for the 60 Ca isotope, 
compared to both Q g -systematics and AA calculations 
(4(±1) x 10~ 13 mb), that makes the search for 60 Ca more 
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TABLE II. Calculated ' Ca production cross section in 
76 Ge, S2 Se + 9 Be reactions. Masses calculated with the 
GXPF1B5 interaction have been used for the Qg-systematics 
and the LISE ++ Abrasion-Ablation model calculations. 



Primary 


Estimation 


59 Ca 


60 Ca 


beam 


method 


cross section 


cross section 






(mb) 


(mb) 


82 Se 


AA 


1.39e-12 


4.83e-13 


82 Se 


Q 3 


1.62e-12 


3.19e-13 


76 Ge 


Q 3 


1.36e-12 


2.47e-13 


82 Se 


EPAX 3 


3.73e-13 


1.65e-14 


76 Ge 


EPAX 3 


3.68e-13 


1.71e-14 



likely to be successful in the near future. 

V. SUMMARY 

The present study of fragmentation of a 82 Se beam 
at 139 MeV/u found evidence for four previously un- 
observed neutron-rich isotopes ( 64 Ti, 67 V, 69 Cr, 72 Mn). 
The longitudinal momentum distributions and cross sec- 
tions for a large number of neutron-rich nuclei produced 
by the 82 Se beam were measured by varying the tar- 
get thickness in a two-stage fragment separator using 
a narrow momentum selection. The momentum distri- 
butions of 126 neutron-rich isotopes of elements with 
11 < Z < 32 were compared to models that describe the 
shape and centroid of momentum distributions. From 
these measurements we have obtained a new set of pa- 
rameters for the semiempirical momentum distribution 
models [Hill. 

The most neutron-rich nuclei of elements with Z = 19 
to 21 have been produced with an enhanced rate com- 
pared to the systematics of the production cross sections 



from the Q g systematics. This trend was previously re- 
ported for the fragmentation of a 76 Ge beam and 
therefore the current results confirm those of our previ- 
ous experiment. This is an indication of a change in the 
nuclear mass surface near Z = 20 for very neutron-rich 
nuclei. This result has been explained with a shell model 
that predicts a subshell closure at N = 34 and a more 
pronounced one at Z = 20. We have shown that pro- 
duction cross section systematics can provide a hint of 
a change of the nuclear mass surface close to the drip- 
line. Plotting the cross section of elements against Q g 
can exhibit sudden changes in slope that are correlated 
with regions of changes in the nuclear structure, such as 
(sub)shcll closures. A correlation to the nuclear mass 
models has been shown via plots of the two-neutron sep- 
aration energy. 

It has been shown that the Abrasion- Ablation model is 
very sensitive to the input mass values for the most exotic 
nuclei. The measured cross sections were best reproduced 
by using masses derived from the full pf shell model space 
with the GXPF1B5 [26| effective interaction modified to 
a recent 54 Ca E x {2^) measurement (28j . 

The cross section for production of 60 Ca using a 
82 Se beam on beryllium has been estimated at 4(±1) x 
10" 16 barn. This estimate is based on both Q g - 
systematics and Abrasion- Ablation calculations using 
masses derived from the shell-model effective interaction 
GXPF1B5 [U]. 
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